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Abstract 

In a previous paper, we described new analytic formulae for optically-thick supercritical accretion flows 
(Watarai 2006, hereafter paper 1). Here we present analytic formulae for optically-thin one-temperature 
accretion flows including the advection-dominatcd regime, using the "semi-iterative" method described in 
paper 1. Our analytic formulae have two real solutions. The first solution corresponds to the advection- 
dominated accretion flow (ADAF), and the second solution corresponds to the radiation-dominated ac- 
cretion flow described by Shapiro, Lightman, & Eardley (the so-called SLE model). Both solutions are 
given by a cubic equation for the advection parameter /, which is the ratio of the advection cooling rate 
Qadv to the viscous heating rate Q V i S , i.e., / = Q a dv/Qvis- Most previous studies assume that / is constant 
(/ ~ 1 for the ADAF). However, it is clear that / should be a function of the physical parameters of the 
radiative-cooling dominated regime. We found that the ratio / can be written as a function of the radius, 
mass accretion rate, and viscous parameter a. Using this formula, we can estimate the transition radius 
from the inner optically-thin ADAF to the outer optically-thick standard disk, which can be measured 
using observations of the quiescent state in black hole X-ray binaries. 
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1. Introduction 

Optically-thin accretion flows have been widely stud- 
ied as a likely origin of the power-law spectral compo- 
nent (up to ~ 100 keV) in Galactic black-hole candidates 
(GBHCs) and nearby low-luminosity active galactic nu- 
clei (LLAGN), including Sgr A* in our Galactic Center. 
Thorne & Price (1975) proposed an optically thin accre- 
tion model for the first time, and Shapiro, Lightman, & 
Eardley (1976; hereafter SLE) modified and applied the 
model to X-ray observations of Cygnus X-l. Their mod- 
els were able to explain the observational properties of 
the low/hard state in GBHCs very well, but the model is 
thermally unstable; i.e. the viscous heating rate cannot 
balance the radiative cooling rate so that the SLE model 
is unstable to small-scale perturbations (see Piran 1978; 
Kato et al. 1998). 

The optically thin advection-dominated accretion flow 
(ADAF) model has been investigated extensively since the 
end of the 1990s (Ichimaru 1977; Narayan & Yi 1994, 
1995a, 1995b; Abramowicz et al. 1995). This model has 
been used to interpret the spectra of black hole X-ray bi- 
naries in their quiescent or low/hard state as an alterna- 
tive to the SLE solution. Since ADAFs have large radial 
velocities, and the infalling matter carries the thermal en- 
ergy into the black hole, advective energy transport can 
stabilize the thermal instability; thus ADAF models have 
been widely applied to explain observations of Galactic 
black-hole candidates (e.g., Narayan et al. 1996; Hameury 
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et al. 1997), the spectral transition of Cyg X-l (Esin et 
al. 1998), and the multiwavelength spectral properties 
of Sgr A* (Narayan et al. 1995; Manmoto et al. 1997; 
Narayan et al. 1997). In addition, many ADAF- like mod- 
els have been proposed including outflows, or convection, 
etc. Since about 2000, these models have been united into 
one category, and began to be called radiativcly inefficient 
accretion flows (RIAFs) (Igumcnshchcv & Abramowicz 
2000; Di Matteo et al. 2000; Abramowicz et al. 2002; 
Narayan 2002; Igumenshchev et al. 2003; Yuan et al. 
2003). 

The self-similar treatment is a powerful tool for inves- 
tigating the flow properties under various physical condi- 
tions without complex numerical computations, for exam- 
ple, the optically thin advection-dominated regime (Spruit 
et a. 1987; Narayan & Yi 1994, 1995), convection- 
dominated flows (Narayan et al. 2000; Abramowicz et 
al. 2002), the optically thick advection-dominated regime 
(Wang & Zhou 1999; Watarai & Fukue 1999; Fukue 2000), 
self-gravity (Mineshige & Umemura 1997), the inclusion of 
magnetic fields (Shadmehri & Khajenabi 2005; Akizuki & 
Fukue 2005), etc. In particular, the self-similar solution 
(hereafter SSS) for the ADAF shows a good agreement 
with numerical solutions except at the inner and outer 
boundaries (Narayan et al. 1997). Thus, the self-similar 
treatment seems to be useful as 0-th order approxima- 
tions. 

To increase the level of approximation in the analyt- 
ical formulae, one more condition is required. One of 
the main problems with the self-similar treatment is the 
energy equation. The energy equation in the canonical 
self-similar solution does not explicitly include radiative 
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cooling, due to the parameterization of / = Q a dv/Qvis, 
where Q a dv and Q v is are the advective cooling rate and the 
viscous heating rate, respectively. Almost all self-similar 
ADAF solutions assume that /=const, so that the radia- 
tive cooling is expressed by Q ra d = (1 — f)Q v is- However, 
/ should be a function of radius, and other physical pa- 
rameters when the radiative cooling is comparable to the 
advective cooling. Moreover if / is close to unity, the ef- 
fect of radiative cooling does not appear in the canonical 
self-similar solutions. Recently, Watarai (2006, hereafter 
paper 1) extended the self-similar solution to the / ^const. 
regime by including the effect of radiative cooling via / 
in the optically thick advection-dominated solutions us- 
ing a "semi-iterative" method, and the derived analytic 
formulae agreed well with numerical solutions over a wide 
range in radius. According to his results, / is given by the 
solution of a quadratic equation, and is a function of the 
mass accretion rate and radius. This raises the question 
of whether it is possible to find a variable-/ solution in the 
optically-thin case using the same method as in paper 1. 
If we can obtain the variable / solution without complex 
numerical procedures, it will facilitate comparison with 
observations of the low/hard state of Galactic black- hole 
candidates or low-luminosity AGNs. 

In this paper, we derive the analytic form of / in 
optically-thin accretion flows and discuss the basic prop- 
erties of the analytic solutions. In the next section, we in- 
troduce the basic equations and the canonical self-similar 
solutions. Our new analytic treatment is presented in sec- 
tion 3. In section 4, we comment on the transition zone 
in bimodal accretion flows. The final section summarizes 
our results. 

2. Self-Similar and Analytic Solutions for 
Optically Thin Accretion Flows 

We introduce the ratio of the advective cooling rate, 
Q a dv, to the viscous heating rate, Q V is, / (= Qadv/Qvis), 
defined in the same way as in previous papers. As for 
an optically thin ADAF, the time scale of energy trans- 
port via Coulomb interaction between ions and electrons 
is longer than that of accreting matter. Because the den- 
sity of accreting gas is very low. Thus, such a flow was 
called an advection-dominated accretion flow, i.e., / ~ 1, 
and it becomes a radiatively inefficient throughout the 
whole of the disk. In this paper we assume that / is not 
constant, but variable. First, the canonical self-similar so- 
lutions are presented in subsection 2.1. Second, using the 
self-similar solutions, we improve the parameterization of 
/ via a semi-iterative procedure described in subsection 
2.3. 

2.1. Basic Equations 

In this subsection, we formulate the basic equations, 
which are the same as those given in paper 1 except for 
the radiative cooling term (see also chapter 8 of Kato et 
al. 1998). We assume a steady axisymmetric accretion 
flow (d/dt = d/dip = 0) and a gas pressure-dominated 
regime. We solve the following hydrodynamic equations 



M = -2-KrT,v r 



in cylindrical coordinates (r, ip, z): mass conservation, 
the momentum equation in the radial direction, angular 
momentum conservation, hydrostatic balance, and the en- 
ergy equation, as follows; 

(1) 
(2) 
(3) 
(4) 
(5) 



dv r 1 dH 



dr £ dr 



dr 



M(£-e in ) = -2irr 2 T rtf 

n 9 9 

Qtis = Qadv + Qr 



ad' 



where M, 



E, £, £- in , and H are the mass accretion 



rate, radial velocity, surface density , £ f=a 2pH, specific 
angular momentum (£~r 2 fl), angular momentum at inner 
edge of the disk, and the scale height, respectively. We 
adopt the Keplarian angular frequency in a Newtonian 
potential, = (GM/r 3 ) 1 / 2 . The r-ip component of the 
viscous stress tensor is expressed by T rip , which is related 
to the total pressure as T rip = —all, where II is the height- 
integrated pressure by II « 2pH , and a is the viscosity 
parameter (Shakura & Sunyaev, 1973). The last term in 
equation (2) represents the correction term for the vertical 
gravity (Matsumoto et al. 1984). We assume we are in 
the gas pressure-dominated regime; thus the equation of 
state is 

P = ~^-pT, (6) 

/!TO H 

where &b, p, and toh are the Boltzman constant, mean 
molecular weight, and Hydrogen mass, respectively. 

Viscosity generates thermal energy through the differ- 
ential rotation of turbulently moving gas, and the viscous 
heating rate is expressed by 

dfl 



Q! 



rT r . 



dr 



(7) 



where To is the linear approximation form of angular ve- 
locity (r$i = — d\n£l/d\nr). The "net" internal energy 
flux at each radius, i.e., advective cooling, is represented 
by the following formula: 



n- M ds M U 



2irr u dr 2nr 2 E s ^ 

The variables To and so represent the temperature and the 
entropy on the equatorial plane, and £ is a dimensionlcss 
quantity, 



1 / dlnpp _ dlapo 
\ d\nr 



r 3 -i 



d\m 



(9) 



where T3 is the ratio of the generalized specific heats (e.g., 
Chandrasckhar 1967), 



r, - 1 = 



/3+i2( 7 -ixi-/?r 

4-3/3 



(11) 
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where 7 is the ratio of specific heats and j3 is the gas- 
pressure normalized by the total pressure. 

By assuming we are in the gas pressure-dominated 
regime (i.e., (3 — > 1) and that po oc r -5 / 2 , and p oc r~ 3 / 2 , 
from the self-similar solutions, then Tx and T 3 tend to 7. 
Finally in gas-pressure dominant regime, £ is given by 



3/ 5/3-7 
? 2 I 7-I 



(12) 



Here the range of 7 is 1 < 7 < 5/3. If 7 is exactly equal to 
5/3, then the flow is iscntropic, and no entropy gradient, 
i.e., Q~ dv = 0. In this paper we set 7 = 1.5; thus £ = 0.5. 

The disk is assumed to be optically thin, gas pressure- 
dominated, with radiation produced only by the free-free 
process: 

Q r "ad = ebremsP 2 ! 11 / 2 ^. (13) 

where £t>rems = 1-24 x 10 21 is the coefficient of the 
brcmsstrahlung emission. The thermal synchrotron and 
its inverse Compton cooling are important cooling pro- 
cesses as in shown by several authors (Kusunose & 
Takahara 1989; Narayan & Yi 1995; Manmoto et al. 
1997; Mahadevan 1997). Such cooling processes are sig- 
nificant near the disk inner region (inside 10 r g ), but 
rapidly decreases at disk outer region (roughly F v oc r~ 2 ). 
Hence we ignore the effect of thermal synchrotron and 
Comptonization in this paper. 

2.2. Canonical Self- Similar Solutions 

To derive the characteristic radial dependence of the 
above equations, we transform the equations (2), (3), (7), 
(8) as follows: 



(14) 



where T v = dlnv r / dhir and Tn = rflnn/dlnr, which can 
be determined from the self-similar solutions. The mo- 
mentum equation in the azimuthal direction is 



Am n K B = 27roTI, 



where B = (1 — l\ n /t) is the boundary term. We also as- 
sume fi = fio^K, where is constant. The energy bal- 
ance equation is given by 



' ad 



v fQvis 



(16) 

where / is constant in canonical treatments. 

The above equations can be solved analytically, which 
leads to the following solutions. We note that all so- 
lutions arc normalized by the critical accretion rate, 
M cr ;t = Le/c 2 = 2ircr g / n es , and the Schwarzschild radius, 
7' g = 2GM/c 2 , where L-^, G, M, c are the Eddington lumi- 
nosity, gravitational constant, black hole mass, and speed 
of light, respectively. The normalized mass accretion rate 
and radius are expressed as to = M / M cr j t and f = r/r g . 



S = S / _1 a _1 mr" 1/2 , 



(17) 
(18) 
(19) 



H 



H Q f 1/2 



T = T ff-\ 

where the coefficient with subscript 
physical quantities: 



"0' : 



(20) 
(21) 
(22) 

represents the 



V rO 



n„ 



H 



1 



(T n - 



3/2)^/ 



-1/2 



V2( 2£ 



\Torio 
Bn c 2 




(23) 
(24) 
(25) 
(26) 

(27) 

(28) 



From the radial dependence of the self-similar solutions, 
we can determine the values of Tn, r„, Tq; i.e., in the 
present parameter sets these constants are given by = 
1.5, £ = 0.5, T v = —0.5, and Tn = — 1.5, respectively. £ and 
fio in equation (12) are then automatically determined, 
and finally all the coefficients can also be evaluated using 
equations (23)- (28). 

Using these self-similar formulae, the advective cooling 
and radiative cooling terms are given by 



<3adv = Qadv/^r 



-l f -3 



rad 

where 



Q?ad«- 2 / 



-2-2 
TO r, 



-lf-B/3 



Ql 



adv 



( 15 ) and 



Ql. 



c 

Kef. 
£br 



^ Br n n 2 



,V2c 



ad 



fc B y 1/2 / 2£ 



(29) 
(30) 

(31) 
(32) 



The parameter dependences of these solutions arc iden- 
tical to those of the canonical self-similar solutions by 
Narayan & Yi (1994) or Wang & Zhou (1999) under the 
condition in which / ^constant. 

2.3. Analytic Solutions under the First Iteration 

In the previous self-similar solutions, the effect of the 
radiative cooling is given by (1 — f)Q v i s - However, the en- 
ergy equation does not explicitly include radiative cooling 
because of the constant / assumption. Using the same 
formulation as in paper 1, we can explicitly include the 
radiative cooling term in /. From the definition of /, 



/ 



' adv 

? + 

v vis 



' adv 



l 



(33) 



' adv 



? rad 
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where g is the ratio between radiative cooling and advec- 
tive cooling. From equation (29)-(32), the explicit form of 
9 is 

9 = 7p~ - go - 12 ^ , (34) 

ad 

where go 



a 2 f 



is the ratio 
-3o-3 



Of QradO tO Q a dvO, i-C, #0 = 

/rado/QadvO oc r^f^g -11 . In this paper the value of 
.9o is go — 1-24 x 10~ 2 for the present parameter sets. 
Substituting equation (34) into equation (33), we finally 
obtain a cubic equation for / in nondimcnsional form: 



f-f + go- 



;l/2 



0. 



(35) 



In the optically thick case, / is independent of a, whereas 
in the optically thin case, / depends on a. This arises 
from the difference in the radiative cooling processes be- 
tween the optically thin and the optically thick regimes. 
The thermal brcmsstrahlung cooling depends on the gas 
density and temperature, so changes in density are implic- 
itly included in / via the viscosity parameter a. 

Consider a cubic equation, x 3 + a\x 2 + a-ix + = 0, 
with real coefficients, a\, 02, and 03. First, we evaluate 
the determinant of the cubic equation, D = Q 3 — i? 2 , as 

-3a2 7-) 2a\— 9aia2+27as 



follows: Q ; 



'-,R- 



54 



The coefficients 

" 2 mf V 2 . 



in equation (35) are ai = —1, a-i = 0, and 03 = g^a 
respectively. When Q and R are real and D > 0, the cubic 
equation has three real roots (Press et al. 1986). In the 
present study, Q — 1/9 and R = (27a3 — 2)/54. Here we 
only consider physically plausible solutions, i.e., < / < 1, 
which results in two real solutions: 



ADAF 



(0) 



1 — 2 cos 



-2tt 



/sle(0) 



1 — 2 cos 



3 

(9 + 47T 



(36) 



(37) 



where 8 is a function of f, m, and a. Here /adaf and 
/sle correspond to advection-dominated and radiation- 
dominated solutions. The explicit form of 9 is given by 



ir.m.a) 



27 mr 1 / 2 
-^-90—3- 



- 1 



(38) 



If we modify / in equations (17)-(22) to /(#), new sets of 
analytic solutions will be obtained. This method is similar 
to an iteration procedure: using the SSS as a first approx- 
imation, we improve the solutions for / via the energy 
equation, which includes the radiative cooling term. We 
note that our present solutions have been obtained by the 
first iteration. Original SSSs (f — 1 for ADAF) in equation 
(17)-(22) are modified by the first iteration (see figure 1, 
2). Further iteration can expect the improvement of the 
solutions but the first improvement will be most signifi- 
cant. In other words, our solutions are still quite rough 
approximation but there is an advantage that is possible 
to describe in the analytical form. 
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Fig. 1. Radial profiles of / for various values of the viscous 
parameter a. Solid and dotted lines represent ADAF-like and 
SLE-likc solutions, respectively. Input parameters are a = 0.1, 
0.3, 0.5, 0.7, and 0.9 from left to right, while the accretion rate 
m is fixed at 0.01. 
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Fig. 2. Same as figure 1, but for m=0.01, 0.03, 0.1, 0.3, and 
1 from right to left. The viscosity parameter a is fixed at 0.5. 



3. New Analytic Solutions 

We derived the analytic form of / considering the ef- 
fect of radiative cooling. Figure 1 shows the solution of 
/adaf((9) and /sle(^) for various values of a. As the val- 
ues of a increase, the optically thin solutions can extend 
to the disk's outer regions. The / value is close to unity 
at large radii for high a. When a is large, the viscosity 
effectively extracts the angular momentum of the matter; 
thus advective energy transport will be important over a 
wide range in radius. However, in the outer regions of the 
disk, radiative losses via brcmsstrahlung emission rapidly 
cool the disk as a decreases, so that ADAF-like solutions 
disappear in the disk's outer region. 

Figure 2 represents the radial dependence of the mass 
accretion rate on /. The optically thin solutions depend 
on the mass accretion rate as well as a. The mass ac- 



Analytical Formulae for Optically Thin Accretion Flows 



5 



cretion rate generally determines the density, which in- 
fluences the cooling rate via thermal bremsstrahlung, so 
that if the density is high, the disk will cool more effec- 
tively. Therefore , the two parameters a and m control 
the basic properties of these solutions through advective 
and radiative cooling. 

Note that figures 1 and 2 have a logarithmic scale on 
the horizontal axis and a linear scale on the vertical axis. 
Therefore, the optically thin ADAF solutions decrease 
rapidly at the disk's outer radii. 

Our results do not conflict with the previous / ~ 1 result 
for low m and high a (Nakamura et al. 1997). However, 
according to Nakamura et al. (1997), ADAF solutions for 
low m and small a, (m, a)=(10~ 5 , 10 -3 ), and for high 
m and large a, (m, a)=(0.1, 0.1), also show / ~ 1 over 
a wide range of radii. This could be due to the outer 
boundary conditions for the ADAF. To discuss the outer 
region of the optically thin accretion flow, it is necessary to 
account for the outer boundary conditions. In some cases 
outer boundary conditions can affect the disk structure 
(Yuan 1999; Yuan et al. 2000); we comment on the outer 
boundary conditions in the next section. 

4. Discussion 

4-1. Transition Radius in Bimodal Accretion Flows 

Some observations of GBHCs indicate an hot, hard X- 
ray component and a cold, soft component are coexistent 
in accretion flow (Thorne & Price 1975, Ichimaru 1977 
for Cyg X-l; Shapiro & Teukolsky 1983; Narayan et al. 
1996 for soft X-ray transients; Esin et al. 1997 for Nova 
MUSCA 1991; Narayan et al. 1997 for V404 Cygni). To 
date, two mechanisms have been proposed to explain this 
hot component: evaporation in the vertical direction of 
the disk (Meyer & Mayer-Hofmeister 1994) or 2-zone ra- 
dial geometry which the outer parts of the disk follow the 
optically thick, while the inner disk switches to the opti- 
cally thin, hot solution, which is secularly stable (SLE). 
The 2-zone radial geometry in accretion flow was solved 
self-consistently for the first time by Honma (1996), who 
pointed out the importance of thermal conduction for a 
smooth transition. With cither mechanism, when radia- 
tive cooling works efficiently, the hot flow shrinks toward 
the equatorial plane, and the optically thin accretion flow 
transitions to the standard disk. Therefore, the energy 
transport in the radial direction is important for a smooth 
transition between the optically thin and thick regimes, 
but it seems that the cooling processes determine the lo- 
cation of the transition radius; this is an important point. 

The location of the transition region between an inner 
optically thin regime and an outer optically thick regime 
can be evaluated using equation (36). If we set /adaf(#) = 
0.5 as a critical value of advection-dominated accretion, 
then equation (36) has a physical solution (0 < / < 1). 
Then, the transition radius is 

The parameter dependence of f tr is consistent with the re- 



sults of Honma (1996), and is very sensitive to the viscous 
parameter a. We have used an independent approach to 
derive this result. Note that equation (39) represents the 
radius where the solution disappears, but Honma (1996) 
obtains the real transition radius. 

Turbulent thermal conduction in a bimodal accretion 
flow gives a finite size for the transition zone (Manmoto & 
Kato 2000a, b). In the transition zone, the heating (cool- 
ing) via turbulent thermal conduction is balanced by ra- 
diative cooling. Turolla & Dullcmond (2000) discussed 
the flow properties at the evaporation radius in the out- 
flow solution by Blandford & Begelman (1999); however, 
they treated the transition radius as an input parameter. 

Recently, Lu et al. (2004) found bimodal solutions 
for relatively high values of a (>0.4) without conduc- 
tion effects. They performed numerical integration from 
the sonic point, and then obtained bimodal solutions. 
However, their method uses a variety of initial conditions 
at the sonic point so that the determination of the tran- 
sition radius is difficult. Here we note that our procedure 
for evaluating the transition radius is purely mathemati- 
cal, and the derived transition radius does not depend on 
the outer boundary conditions. 

Figure 3 represents the radial structure of a bimodal 
accretion flow based on our model. The inner region is 
described by the optically thin solution, and the outer re- 
gion is a solution with which the standard optically thick, 
geometrically thin accretion disk by Shakura & Sunyaev 
(1973) is connected by the same accretion rates and a pa- 
rameters as the inner region. The inner boundary term 
in the standard model, 1 — \/r- m /r, was disregarded. Here 
the ri n represents the disk inner radius. We ignore radial 
energy transport between the optically thin regime and 
the optically thick regime. In our solutions, the transition 
radius is not an input parameter unlike previous studies. 
Our estimate of the transition radius is very sensitive to 
the angular frequency, with f tr cx g^ 2 ex fig, and so we can- 
not compare our results with previous papers at this stage. 
Our model is very simple, so it will be useful for discussing 
the rough structure of the flow; however, numerical cal- 
culations will be required for detailed comparisons with 
observations. This will be carried out in future work. 

4-2. Effect of Outer Boundary Conditions 

In deriving a solution for optically thin accretion flows, 
it is necessary to handle the initial value problem with care 
because the outer boundary conditions can sometimes 
change the whole structure of the accretion flow. For ex- 
ample, a Kcplcrian angular frequency at the outer bound- 
ary was assumed in several early studies (Narayan et al. 
1997) but Nakamura et al. (1996, 1997) assumed a sub- 
Keplerian angular momentum and a particular tempera- 
ture at the outer boundary (see also Chen et al. 1997). In 
extreme cases, the outer boundary condition dramatically 
changes the entire disk structure, for example, when the 
outer boundary condition is set to high mach number and 
sub-Kcplcrian angular velocity (Yuan 1999; Yuan et al. 
2000). For this reason, it is not meaningful to compare 
the new analytical solutions with numerical solutions. 
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Fig. 3. Radial dependence of bimodal accretion flows with 
various a. Panels arc the surface density (top left), radial 
velocity (top right), gas temperature (bottom left), and scale 
height (bottom right), respectively. We adopted the optically 
thick standard disk solutions in the disk outer region. 



We suggest that the plausible outer boundary condi- 
tions must be determined from observations of the disk's 
outer regions in the objects being studied. 

4-3. Stability of Transition Zone 

The time-dependent evolution of the transition zone, 
including the effects of turbulent energy transport, was 
recently computed by Gracia et al. (2003). They con- 
firmed the result of Manmoto et al. (2000a), Manmoto 
& Kato (2000b), and found that the transition zone is 
highly variable, suggesting that it may be the origin of 
the quasi-periodic oscillations observed in Galactic X-ray 
binaries. However, if the magnetic field in the transition 
zone is taken into account, the gas density is confined by 
the magnetic field, which can cause the cooling energy at 
a given radius to change and potentially alter the location 
of the transition radius. 

Changes in the transition surface due to magnetorota- 
tional instability (MRI) were examined by Kato (2000) 
using linear analysis. Linear MRI perturbations gener- 
ated in the outer standard disk penetrate into the in- 
ner ADAF, which results in a highly variable transition 
zone. Recently, analytic solutions, including magnetic 
fields, have been investigated extensively by several au- 
thors (Meier 2005; Shadmchri & Khajenabi 2005; Akizuki 
& Fukue 2006). According to Machida et al. (2006), the 
magnetically dominant configuration survives for much 
longer than the thermal time, and is comparable to the ac- 
cretion time scale, but it is shorter than the viscous time. 
However, such a long time calculation is still rare, so the 
result has not been reached consensus yet. We hope to be 
confirmed for a long time by the simulation a lot near the 
future. 

It is possible that oscillation of the transition surface 
produces the quasi-periodic oscillations commonly ob- 
served in Galactic black hole candidates (see also Kato 
2000). This situation is similar to oscillations of the 
Alfvcn radius of the magnetosphere around a neutron star 
(Shapiro & Teukolsky 1983). Thus MHD simulations of 
the region around the transition surface will be an intrigu- 
ing topic for further research. 

5. Conclusion 

We described analytic formulae for accretion flows in 
the optically thin regime. The ratio of the advection cool- 
ing rate to the viscous heating rate is given by the solu- 
tion of a cubic equation, and is a function of f, m, and a. 
Solving the cubic equation, we can obtain not only ADAF 
solutions, but also SLE-likc solutions in a simple way. Our 
new solutions are physically valid, because both solutions 
originate from same cooling processes. The transition ra- 
dius can be evaluated, where f{9) = 0, and is consistent 
with Honma's (1996) result. We expect that our analytic 
technique will be useful for order estimations of optically 
thin accretion flows; moreover, it will be useful to provide 
initial conditions for numerical simulations. 
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